Influence of flow confinement on the drag force on a static cylinder 
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The influence of confinement on the drag force J" on a static cylinder in a viscous flow inside a 
rectangular slit of aperture ho has been investigated from experimental measurements and numer- 
ical simulations. At low enough Reynolds numbers, F varies linearly with the mean velocity and 
the viscosity, allowing for the precise determination of drag coefficients A|| and A_l corresponding 
respectively to a mean flow parallel and perpendicular to the cylinder length L. In the parallel 
configuration, the variation of A[| with the normalized diameter /3 = d/ho of the cylinder is close 
to that for a 2D flow invariant in the direction of the cylinder axis and does not diverge when 
(3 — 1. The variation of A|| with the distance from the midplane of the model reflects the parabolic 
Poiseuille profile between the plates for /3 ^ 1 while it remains almost constant for /3 ~ 1. In the 
perpendicular configuration, the value of Ax is close to that corresponding to a 2D system only if 
/3 <C 1 and/or if the clearance between the ends of the cylinder and the side walls is very small: in 
that latter case, Ax diverges as /3 — > 1 due to the blockage of the flow. In other cases, the side flow 
between the ends of the cylinder and the side walls plays an important part to reduce Ax: a full 3D 
description of the flow is needed to account for these effects. 

PACS numbers: 47.15.G-, 47.15. Rq, 47.60.Dx,47.63.mf,47.11.Fg 



I. INTRODUCTION 

Because of its fundamental and practical implications, 
considerable research efforts have been devoted to the 
study of fluid flow past fixed or moving slender bodies. 
Typical examples include the settling of suspensions of 
solid particles such as found in the paper industry or the 
injection of fibers P, [^. Existing literature reports that, 
for flows in confined geometries such as pipes, hydrody- 
namic interactions between the fibers and the wall have 
a great influence on the fiber orientation and, in turn, on 
the flow properties Q. Before addressing such complex 
situations involving a large number of particles, simple 
models considering isolated fibers have recently been de- 
veloped 0]. In addition, the recent development of appli- 
cations involving micro and nano rod-like objects p, Q 
also raise questions on the hydrodynamic forces on ob- 
jects placed in confined geometries such as microfiuidic 
channels. 

In this paper, the hydrodynamic forces F acting on a 
static cylindrical rod inside a viscous flow in a slit of rect- 
angular cross section ho x W (with ho -C W) are deter- 
mined both experimentally and numerically. We have, in 
particular, compared the cases of a cylinder parallel and 
perpendicular to the flow. The rod is a cylinder of high 
aspect ratio, i.e. its length L is always larger than its 
diameter d {L ^ d) but its length is of the same order as 
the slit size W. The effect of the confinement due to the 
two closest plane plates of the slit is particularly investi- 
gated: more precisely, the influence of the ratio f3 — d/ho 
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is studied over a broad range of values up to /3 '--^ 1 (very 
strong conflnement). The influence of the distance W 
between the two lateral sides is also analyzed: it is par- 
ticularly significant for cylinders normal to the mean flow 
and blocking it partly, resulting in large hydrodynamic 
forces. 

For viscous flows (either confined or not), the drag 
force F per unit length is proportional to the dynamic 
viscosity r/ and to the velocity U far from the cylinder 0] ■ 
For unbounded flows, the leading term of the proportion- 
ality coefficients, named in the following A|| = F\J{rjLU) 
(resp. A_L = F±/{r]LU)) for flow respectively parallel 
and perpendicular to the axis of the cylinder, is of order 
e=l/\n{d/L) 

For slightly confined flows, wall correction terms in- 
crease the drag force. The configuration in which a small 
cylinder sediments half-way between parallel vertical 
plates separated by a large distance ho (d ^ L <C /iq < 
W) has been studied experimentally by de Mestre [IQ] . 
In this case, one has A = ai{e + (a2^ + as)^^), where 
the parameters ai are constants, depending only on the 
orientation of the cylinder which is either vertical or hor- 
izontal. This configuration has also been studied in the 
limit of cylinders of length large compared to the aper- 
ture {ho <C i <C W). In this case, the influence of the 
boundaries on the drag is dominant and A depends solely 
on the ratio /3 = d/ho (and not on the cylinder length L) 
and scales like l/ln(/3) for /3 < 1 P, [II, [H . 

These theoretical predictions were confirmed and ex- 
tended to large values of (3 by 2D numerical simula- 
tions [14J and experimental measurements [l^, [l^, 
in the case of cylinders moving with their axis normal to 
the flow. The situation where the cylinders are flxed 
have been also considered in the limit of low J9,, 
and high [m, |F 
complex fluids 
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ever scarce in this geometry: Dhahir and co authors (26[ 
measured forces on a cylinder of low aspect ratio (d ~ L) 
for different fluid rheologies while Rehimi et al [23] use a 
geometry similar to ours, but do not perform force mea- 
surements. 

In the present experiments, forces are measured on a 
static cylinder in a long slit of rectangular section where 
the flow takes places. Both /? = d/ho and L/W are 
varied as well as the distance of the axis of the cylinder 
from the center plane of the slit; both cases of a cylinder 
parallel and perpendicular to the mean flow are studied. 
In the perpendicular case particular attention is given to 
the influence of the flow between the ends of the cylinder 
and the sides of the slit (if L < W): flow in this region is 
highly three dimensional and 3D simulations are needed 
to estimate it. 

The other, parallel, configuration has not been studied 
up to now to our knowledge. However, previous authors 
studied theoretically in the viscous regime [28j and exper- 
imentally in the inertial regime [1^ the related problem 
of the forces on a cylinder located inside another, coax- 
ial, one. In the viscous regime, the velocity of a cylin- 
der falling inside another one has also been investigated 
[sot . this parallel case, we measure the drag forces 
induced by the flow on cylinders of different diameters 
(0.04 < (3 = d/ho ^ 0.83) and for different locations in 
the aperture of the slit. In order to extend the range 
of physical parameters investigated, 2D numerical sim- 
ulations are performed; they allow, in addition to dis- 
criminate between the contributions of the pressure and 
viscous shear forces to the global measured drag force. 

The experimental setup is presented in section |TT] and 
the numerical method in section IIIII The experimental 
and numerical results are reported in sections IIVI and |V] 
respectively for cylinders parallel and perpendicular to 
the flow. 



II. EXPERIMENTAL SETUP AND 
PROCEDURE 

The experimental setups used for the determination 
of A|| and are shown in Fig. [T] they consist of a 
slit of rectangular cross section placed vertically and 
made of two transparent milled polymethyl methacry- 
late (PMMA) plates. The cell aperture has the constant 
value ho = 4.9 or 0.75 mm except in the upper 60 mm 
of the cell where it increases with height from 0.75 mm 
(resp. 4.9 mm) to 5 mm (resp. 10 mm) for the first (resp. 
second) model. This Y-shapc profile was designed to ease 
the insertion of the cylinders into the cell. 

The cylinders are hung on the hook of computer con- 
trolled scales (Sartorius CP 225) measuring drag forces 
of the flowing fluid on the cylinder in a range of 10~^ 
to 8 X 10~^N. A precision translation stage allows one 
to displace the cylinders across the gap of the cell: in 
this way, the hydrodynamic forces can be measured as 
a function of the distance from the walls. Furthermore, 




FIG. 1: Left: schematic view of the experimental setup used 
to measure A||, L £ [49, 200] mm, I = 60mm, W = 90mm, 
ho is either 0.75 or 4.9 mm. Right: setup for measuring A^, 
ho = 4.9 mm, L G [20, 89] mm. 



the latter are transparent allowing one to determine pre- 
cisely the location of the cylinders within the cell. For the 
cell with the largest aperture, side views can also be ob- 
tained so that the distance separating the cylinder from 
the walls may also be measured; this also allows one to 
control the parallelism of the object with respect to the 
wall. 

A gear pump sucks the fluid at the bottom side of the 
cell and reinjects it into a bath covering the top side (the 
fluid flows therefore always vertically downwards). The 
fluids are either pure water or water-glycerol mixtures 
with a relative mass concentration of glycerol ranging 
from 50 to 80%. 

The value of the viscosity t] is determined by first mea- 
suring the density p of the solutions and their tempera- 
ture before each series of experiments by means of an 
Anton Paar 35N densimeter. The viscosity rj is then de- 
termined from the measured density and temperature by 
interpolating tabulated values. The final uncertainty on 
the value of the viscosity is about 3% and is lower for all 
other parameters. The relative influence of viscous and 
inertial effects is characterized by the Reynolds number 
Re ~ phoU /rj in which U is the mean velocity far from 
the cylinder. 

Measurement of A|| were performed for several cylin- 
ders with different diameters (see Tab. HI, which were 
either rigid (glass, copper, iron, PMMA) or flexible 
(polyester or silk threads). The flexible threads were 
stretched prior to the experiments in order to remove 
their residual curvature. These threads include multiple 
fibers so that their diameter is not constant and varies 
periodically. Yet, such variations were found to have a 
negligible effect on the drag force: in the next parts, these 
threads are characterized by their mean diameter. 
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TABLE I: Experimental parameters corresponding to the 
measurements of A||. d, L: diameter and length of the 
cylinder; ho: cell aperture in the constant aperture region; 
P = d/ho; T]: dynamic viscosity of the solutions. The symbols 
characterizing the different experiments are the same as in 
the experimental figures. For data corresponding to symbols 
(□), (ffl) and (H), L is the difference between the lengths of 
two cylinders. 





d(mm 


) /t(mm) 


13 


L(mm 


) ?7(mPa.s) 


symbol 


glass 


f.5 


4.9 


0.31 


89 


40.0 


□ 


PMMA 


3.2 


4.9 


0.65 


110 


22.4 


ffl 


PMMA 


4.05 


4.9 


0.83 


100 


21.7 




optical fiber 


0.14 


4.9 


0.029 


195 


35.0 


< 


polyester 


0.20 


4.9 


0.041 


177 


17.5 


O 


polyester 


0.20 


4.9 


0.041 


177 


7.34 


V 


silk 


0.45 


4.9 


0.092 


153 


37.2 


> 


glass 


1.5 


4.9 


0.31 


151 


6.5 


A 


glass 


1.5 


4.9 


0.31 


138 


40.0 


o 


iron 


2.0 


4.9 


0.41 


177 


32.0 


1X1 


iron 


4.0 


4.9 


0.82 


184 


37.6 





PMMA 


4.05 


4.9 


0.83 


179 


21.7 


® 


polyester 


0.18 


0.75 


0.24 


135 


119 


T 


silk 


0.45 


0.75 


0.6 


82 


24.0 


► 
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FIG. 2: Schematic 2D representation of the experimental con- 
figuration, (a): parallel case, (b): perpendicular case 



value of U almost coincide: this demonstrates the very 
good reproducibihty of the measurements and the lack of 
hysteresis between the phases during which the flow rate 
is increased or decreased. 



TABLE IL Values of the experimental parameters corre- 
sponding to the measurement of Ax. All these experiments 
were performed on the cell of aperture ho — 4.9 mm and width 
W = 90 mm; Re < 5. 





d(mm) 


/3 


range of L(mm) 


symbol 


position 


steel 


0.98 


0.20 


85 


T 


center 


glass 


1.5 


0.31 


37.5 - 88 


▲ 


center 


brass 


2.96 


0.60 


86.5 


■ 


center 


PMMA 


3.2 


0.65 


44.7 - 88 


♦ 


center 


PMMA 


4.05 


0.83 


67.5 - 88.7 


• 


center 


glass 


1.5 


0.31 


37.5 - 84 


A 


wall 


brass 


2.96 


0.60 


86.5 


□ 


wall 


glass 


4.15 


0.85 


25.4 - 88.8 


o 


wall 


glass 


4.15 


0.85 


57.6 


+ 


wall 



III. NUMERICAL SIMULATION PROCEDURE 

In order to determine numerically A|| and A^, the 
Stokes and incompressibility equations must be solved 
with appropriate boundary conditions. The shear stress 
and the pressure at the surface of the cylinder are then 
computed and added in order to determine the total hy- 
drodynamic force. For flow parallel to the cylinders, 2D 
simulations provide reliable results; in the perpendicular 
configuration, 2D simulations have also been used but 
3D simulations are more realistic in several cases. Fig- 
ures [5] displays the geometry used in the 2D simulations: 
the pressure gradient is applied in the z direction, which 
is along the cylinder axis for A|| (see Fig. and per- 
pendicular of the cylinder axis for Aj^ (see Fig. ^jp). 



For measuring X± (see Tab. lll| . the rigid cylinders are 
placed horizontally in the cell of largest aperture {i.e. 
ho = 4.9 mm); their center is halfway between the side 
walls of the cell. There are hung using threads of small di- 
ameter (100 /im) as shown in the right drawing of Fig. [1] 
Therefore, the drag forces on the threads and on the ver- 
tical rod add up: the former is however generally small 
compared to the latter. Moreover, this extra force may 
be estimated numerically (see Sec. IIII[) and subtracted 
from the measurements. 

The flow rate is increased stepwise from to Q = up to 
Q = 400mL/min (resp. Q = 1400mL/min) for the water 
glycerol mixtures (resp. for water), and then decreased 
down to Q = 0. Three such cycles are performed in order 
to verify the reproducibility of the measurement. Figures 
[3] and [3] display the variations of the drag force (averaged 
during each constant flow rate step) as a function of the 
mean velocity U. All data points corresponding to a same 



A. Computation of A 

For a constant pressure gradient parallel to the z 
axis, the velocity is everywhere parallel to z and V = 
V{x,y)ez due to the translational symmetry of the sys- 
tem. The governing equation of the flow then reduces to 
a 2D Laplace equation, which has been solved by means 
of the finite element program FreeFem-!--!- [sij . The grids 
contain at least 12000 nodes (values of the force com- 
puted with a finer mesh size were identical to within 2%). 
The hydrodynamic force per unit length / is computed 
by adding to the shear force the contribution of the pres- 
sure {j/3'^dp/dz , where /? = d/ho). The influence of 
the lateral side plates is taken into account by setting a 
zero-velocity boundary condition for x = ztW/2. A cell 
of infinite width is modeled by assuming, for x = ztW/2, 
a parabolic variation of the velocity with the distance y 
from the midplane. 
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The force measured experimentally is the integral of 
the forces acting along the full length of the cylinder. 
Assuming that the velocity profile in the gap of the cell 
depends only on its local aperture h{y) (lubrication ap- 
proximation), the total force is computed numerically 
from: 

FiubXL) = I ^ f{z)dz, (1) 

in which L is the length of the part of the cylinder located 
inside the constant aperture domain (See Fig. [1]) and I 
corresponds to the part inside the Y-shaped section (in 
the upper fluid bath, the fluid velocity is small enough 
so that its contribution to the force can be neglected). 



TABLE III: Comparison between values of X± for a cylinder 
halfway between the plates from the present work and from 
previous numerical and theoretical studies. The analytical 
solution of ref. uses Eq. ([21; all the results of the table 
correspond to 2D systems. 

Numerical Numerical Analytical Analytical 
/3 Present work Richou(2004) Faxen (1946) Richou(2004) 



0.01 


5.100 


5.309 


5.109 




0.1 


13.34 


13.74 


13.36 


11.52 


0.4 


72.69 


73.28 


72.93 


72.55 


0.6 


262.4 


266.8 




265.3 


0.8 


1850 


1884 




1866 


0.96 


1.205 X 10^ 


1.149 X 10^ 




1.208 X 10^ 


0.99 


3.955 X 10® 


3.174 X lO'' 




3.965 X 10"^ 



B. Computation of Ax 

1. 2D computation 

Here, we consider the configuration in which the flow 
is normal to the cylinder (see Fig. [^b). Far from the 
cylinder, the velocity has a parabolic Poiseuille profile. In 
this 2D approach, the flow is assumed to be invariant in 
the X direction so that: V = Vy{y, z)ey + Vz{y, z)ez- The 
2D Stokes equation is solved numerically by means of the 
FreefemH — h program, using a mesh containing at least 
15000 nodes. The length of the computational domain in 
the z direction is 8 times its height in the y direction (we 
checked that choosing a longer computational domain in 
the z direction has a negligible influence on the value of 
the force). The shear and the pressure forces are then 
computed from the velocity field. 

In order to validate the present method, we compared 
its predictions in the particular case of a cylinder located 
halfway between the plates to those available in the lit- 
erature. In particular, Richou et al. ^ estimated ana- 
lytically an asymptotical value of in the limit when 
the free space between the cylinder and the walls is small 
(/? « 1) by using the lubrication approximation. This 
approximation is valid in regions where one of the front 
walls and the surface of the cylinder are close to each 
other and nearly parallel. Since the volume flow rate 
is conserved along the stream-tubes, these are the same 
regions in which the fluid velocity is the highest (and 
therefore also the friction force and the pressure gradi- 
ent). Rewriting these results 4] with our notations gives: 



Ai(/3,5 = 0) 



6ttV2- 



1/2 



(1-/3) 



3/2 



f3' 



3/2 



(2) 



where S is the distance between the middle of the cell and 
the axis of the cylinder (See Fig. Oa). The pressure term 
is dominant when (3 — > 1. Still using the lubrication 



approximation, we extended this result to the case of a 
cylinder touching the wall in the constant aperture region 
(thickness of the free space equal to zero on one of the 
sides) with: 



X±iPJ={ho-d)/2) 



1 



2V2 



Ai(/3,<5 = 0). (3) 



Table Hill compares, for different values of /3 = d/ho, the 
results of the present 2D numerical simulations to the 
predictions of Eq. Q and to the theoretical and numeri- 
cal results of other authors for a cylinder located halfway 
between the plates. The good agreement between the dif- 
ferent values validates the present numerical procedure. 



2. 3D computation 

The 2D approaches described above assume that the 
cylinder is infinitely long: they describe correctely there- 
fore only the case of a cylinder of length equal to the 
width of the slit {L ^ W in Fig. [T]). U L<W, there is a 
deviation of the flow lines from the z = est planes (Fig. [J) 
resulting from the free space separating the edge of the 
cylinder and the lateral walls. If this lateral clearance 
is large, this reduces strongly the force on the cylinder 
compared to the 2D configuration. 

In order to estimate the magnitude of this latter ef- 
fect and compare the results to the experimental obser- 
vations, we solved numerically the SD Stokes equation in 
some of our experimental configurations by means of the 
FreefemSD program (33 |. The numerical discretization 
consists in a Pl-Pl finite element method using penalty 
for the pressure stabilization; the number of vertices is 
at least 70000. The ratio between the numerical lengths 
in the directions x and y is 18 (like in the experimental 
cell), and the ratio between the numerical lengths in the 
directions y and z directions is 12. 
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IV. EXPERIMENTAL AND NUMERICAL 
VARIATIONS OF A I 

Most previous experimental studies investigated only 
configurations in which the cylinders are weakly confined 
{i.e. (3^1) and/or located halfway between the walls of 
the channel. The present experimental setup has allowed 
us to investigate the cases in which the confinement is 
strong as well as the variations of the force when the 
cylinder comes close to the wall. 



A. Variation of the drag force on a cylinder 
parallel to the flow with the mean velocity 

In the present study, the dependence of the drag force 
on the geometry of the system for viscous flow is ad- 
dressed. For this flow condition, the drag increases linarly 
with the mean flow velocity, so that A|| is independant of 
the velocity. Therefore, before analyzing the dependence 
of A|| on the geometrical parameters of the flow, we in- 
vestigated the domain in which the force and the velocity 
are proportional for fluids of different viscosities. 



600- 



S 400- 
o 

LL 



200 




5 10 
Velocity (mm.s"^) 



15 



FIG. 3: Variation of the force on a vertical cylinder located 
halfway between the plates as a function of the mean velocity 
(7 of a water glycerol mixture, d = 1.5 mm, ho = 4.9 mm, 
f3 = 0.3 and rj = 40.0mPa.s, Reynolds number Re < 2. 2. (A): 
experiments. Dashed line: numerical computation. 

Figure [3] displays a set of measurements obtained using 
a water-glycerol mixture as the flowing fluid: the dashed 
line is the variation of the force with the velocity pre- 
dicted by the 2D numerical simulations of section IIII Al 
The difference between the experimental data and the nu- 
merical results is lower than 5% for all the experiments 
performed in the parallel case, without any adjustable 
parameter. The linear increase of the force with the fluid 
velocity implies that the viscous effects are dominant. 
This is in agreement with the low value of the Reynolds 
number Re which is less than 3 for all the experiments 
using the water-glycerol solution including that of Fig.[3l 

This condition is not satisfied any more for water: in 
this case, the Reynolds number can be as high as 200 



80 



3 60 

o 
o 

o 

u. 40 



20 



A , 



A, 



10 20 30 

Velocity (mm.s"^ ) 



40 



FIG. 4: Same measurements as in Fig. [3] but using water 
(77 = l.OmPa.s): Re ~ 220 for the highest flow rate. (A): 
experiments. Dashed line: numerical computation. 



and inertial effects are important even at moderate ve- 
locities as can be seen in Fig. |4l The force follows at 
first a linear trend with U : moreover, the slope at the 
origin is in good agreement with the theoretical value 
(dashed line). However, and in contrast with the case 
of the water-glycerol solution, the curve deviates from 
the linear trend at higher velocities which reflects the 
increase of the inertial effects: the deviation becomes 
large for mean velocities U > lOmm.s""'^ corresponding 
to Reynolds numbers Re > 50. In the following, only 
results corresponding to the domain of linear variation of 
F with U are reported. 



B. Variation of A with the diameter 

In this section, we study how confinement influences 
the hydraulic forces on cylinders with their axis vertical 
{i.e. parallel to the flow) and located halfway between 
the two parallel vertical walls of spacing hg. For this 
purpose, we determine A|| from force measurements and 
we study its variation as a function of the ratio (3 = d/ho 
by using cylinders of different diameters d. 

However, in this configuration, the local distance be- 
tween the front walls varies continuously in the Y-shaped 
section at the top of the cell. As a result, the drag force 
per unit length increases continuously along the rod and 
is only constant in the region of constant aperture /iq. In 
order to determine the value of the drag specifically in 
this latter region, we repeat the experiment twice with 
cylinders of two different lengths Li and L2 (large enough 
to reach the constant aperture domain), but otherwise 
identical. From Eq. ([1]) , the difference between the forces 
measured at a same mean flow velocity U on the two 
cylinders is: AF — X^rjU {Li ~ L2). 

Figure [5^ displays (solid lines) the variation of A|| 
with (3 predicted from the 2D numerical simulations, 
for a cell of aspect ratio W/h^ = 18 {i.e. of aper- 
ture ho = 4.9 mm). Values obtained by subtracting 
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0.0 0.2 0.4 0.6 0.8 p 1.0 



FIG. 5; Variation of A|| vs. f3 — d/ho for cylinders located 
halfway between the walls. Symbols of experimental data 
points are defined in Tab. ID Solid lines: 2D numerical pre- 
dictions, (a): cell of aspect ratio W/ho « 18 {ho = 4.9mm). 
Dashed line: numerical value in the limit W/ho oo. (b): 
cells of very large aspect ratio W/ho > 120 {h = 0.75mm); 
dotted line: viscous shear force term, dashed line: pressure 
term. 



the forces measured for two cylinders of different lengths 
((□jffljKl) symbols) are in good agreement with the nu- 
merical model. 

In the present experiments, the part I of the length 
of the cylinder inside the Y-shaped section is generally 
much smaller than the length L inside the constant aper- 
ture zone. In this case, we estimate A|| from a single 
measurement of the force F by assuming that, although 
the local aperture h{z) varies, A|[ is the same in all sec- 
tions and that the corresponding local mean velocity is 
U{z) = ho U/h{z) (in order to insure mass conservation). 
Eq.© becomes then: F = X\\7]U{L + J°^{h / h{z))dz) 
with h{z) = ho + (Hq — hi)z/l {hi is the aperture at 
the top of the Y-shaped section). A|| is then related 
to the force F by: A|| = F/{rjUL*) in which L* = 
L + I ho ln{j^) / {hi — ho) is an equivalent length. All 
data points, except those corresponding to the (□, ffl, Kl) 
symbols in Figure [5^, were obtained by this "equivalent 
length" method: its validity is demonstrated by the small 
difference between values obtained for same experimental 
parameters by the two methods {® and H symbols). 



Points (O) and (v) correspond to the same experimen- 
tal configuration except for the viscosity which differs by 
a factor of 2: they coincide within experimental error 
which confirms that the force is proportional to the vis- 
cosity. 

In figure[5^, experimental data points are (as expected) 
closer to the numerical values taking into account the fi- 
nite aspect ratio W/ho = 18 of the cell (solid line) than 
to those assuming an infinite width W (dashed line). 
The difference between the two curves is however only 
of the order of 10%, indicating that the effect of the lat- 
eral boundaries is weak. This correction becomes com- 
pletely negligible for the narrower cell {ho = 0.75 mm 
and W/ho — 120) and the results are then the same as 
for W/ho ^ cxD. In this case, it is more difficult to control 
experimentally the position of the cylinders: the two ex- 
perimental values obtained (using the equivalent length 
approach) are, however, in good agreement with the nu- 
merical predictions (see Figure [S}d) . 

Still in the high aspect ratio limit, we computed nu- 
merically the pressure force term (dashed line in Fig.[S)3) 
and the shear force term (dotted line in Fig. [SJd). On 
the one hand, for /3 > 0.2, the viscous contribution lev- 
els off and decreases for (3 > 0.8; on the other hand, the 
pressure contribution increases sharply when (3 increases. 
For /3 > 0.2, the sum of the two contributions {i.e. \\\) 
increases almost linearly with (3 (A|| « 2.1 -I- 13.8 /?), and 
does not diverge for (3=1, due to the weak perturbation 
of the flow. 



C. Variation of A with the location of the cylinder 
in the aperture 




-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 



FIG. 6: Variation of A|| with the normalized distance 5 /ho 
from the middle plane {\S/h\ < |0.5 - d/{2ho)\ = 0.35) for 
d — 1.5 mm and ho — 4.9 mm. Other parameters are the same 
as for the (A) symbols in Tab. U] (□, O, o): repeat experi- 
mental runs for the same parameter values as (A). Solid line: 
numerical computation of the force; dashed line: Poiseuille 
profile. Insert: variation of the ratio of the values of A|| at 
the wall and at the center as a function of /3. (A) and (ixi): 
experimental values; solid line: numerical computations. 
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The experimental setup also allows one to move pre- 
cisely the cylinder across the gap of the cell. The nor- 
malized distance S/h^ along y of the axis of the cylinder 
from the middle plane of the cell is then determined with 
a precision better than 50 /im using the side view pic- 
tures. The experimental values of A|| arc plotted as a 
function of 5 /h^ in Fig. [6] : in this case, the sequence of 
displacements towards and away from the walls is per- 
formed three times. For all the values of (3 considered in 
the present work, this curve has a parabolic shape with 
a maximum for 5 /ho = 0. The coefficient A|| decreases 
when approaching the wall (solid line Fig. [5]) but less than 
if it followed a parabolic variation for a Poiseuille veloc- 
ity profile with no cylinder present (dashed line). The 
influence of the finite diameter of the cylinder, particu- 
larly when it is of the order of the distance to the wall, 
accounts for this difference. 

The insert of Figure [S] displays the variation with (3 of 
the ratio of the values of A|| for a cylinder in contact with 
the walls and in the middle of the cell. This curve is well 
fitted by the function 1 — 0.9 x exp(— 3.4/3) in the range 
0.03 < /3 < 1. Even though the overall variation is always 
parabolic, the decreasing trend of A|| near the wall is 
stronger for lower values of /3; for /3 ^ 1, this variation is 
expected to be similar to the Poiseuille parabolic profile. 



V. EXPERIMENTAL AND NUMERICAL 
VARIATIONS OF Ax 

While the cylinder perturbs only weakly the flow when 
its axis is parallel to it, the perturbation is much larger 
in the perpendicular configuration. More precisely, the 
flow section may be signiflcantly reduced in the vicinity 
of the cylinder, particularly as the normalized diameter 
(3 = d/ho 1. This blockage effect forces the fluid to 
flow around the cylinder: it increases the local velocity 
and pressure gradient and, therefore, the drag force F. 
Paragraph IV Al reports measurements of F for different 
values of /3 and for various locations of the cylinder in the 
slit section (but always for L ~ W). As in the parallel 
case (see section IIV Ap , only experiments in which the 
drag force is proportional to the velocity are used (this 
corresponds to a Reynolds number lower than 5). 

The blockage effect discussed above is reduced when 
the cylinders do not span the full width W of the slit {i.e. 
L < W): then, a part of the flow takes place between the 
ends of the cylinder and the side walls. This bypass ef- 
fect reduces in turn the drag force F. Paragraph IV Bl 
reports experimental measurements of the variation of 
Aj_ as a function of L/W together with the 3D numeri- 
cal simulations which are needed to reproduce the highly 
three-dimensional bypass flows. 

We observed that, when perpendicular to the flow, the 
cylinders move towards the middle of the cell while re- 
maining parallel to the plates (except for short cylinders) 
even for Reynolds numbers as low as 0.1. Such a repul- 
sion effect has already been reported by [2l|, [S^l and is 



related to inertial effects of small magnitude (in a pure 
viscous flow, the lift would be zero [3). In contrast, 
the studies of Zovatto and Juarez [2ll.T22| predicted an 
attraction by the walls due to variation of the shear in 
the slit gap, resulting in a shift - increasing with the 
Reynolds number - of the equilibrium position toward 
the walls. However, this phenomenon was observed only 
for Reynolds number larger than 80 while the present 
study deals with Reynolds number {Re < 5); this effect 
may therefore be expected to be negligible in the present 
work and the cylinders should reach (as indeed observed) 
an equilibrium in the middle plane of the cell. 



A. Variation of Ax with the diameter and the 
location of the cylinder in the aperture 




^ 1 1 1 1 r—^ 

0.0 0.2 0.4 0.6 0.8 „ 1.0 



FIG. 7: Variation of Ax as a function of /3 for L ~ W. 
Cylinder in contact with a wall: open symbols (experiments); 
dashed dotted line {2D numerical computations). Cylinder 
in the middle plane of the cell: filled symbols (experiments); 
solid line {2D numerical computations, see Tab.lng. Dotted 
line: ratio of the values corresponding to the solid and dashed 
dotted lines. 

In these experiments, the length of the cylinder is cho- 
sen as close as possible to the width of the cell (0.93 < 
L/W < 1) in order to minimize the lateral bypass flow 
(this point is discussed in detail in the next subsection). 
In this case, Figure [7] displays the variation of Aj_ as a 
function of /?, both for cylinders located in the middle 
plane of the cell and in contact with the front wall; this 
latter case is achieved experimentally by inserting a mag- 
netic wire in the cylinder and attracting it with a magnet. 

Here, too, the experimental results and the 2D numeri- 
cal simulations agree (to better than 10%). Experimental 
values of Aj^ up to 2500 are measured while, as shown in 
Fig. [5l the maximal value of A|| is about 18; such large 
values and the divergence of Aj_ when (3^1 are due to 
the small gap left for the flow between the cylinder and 
the front walls. As this gap decreases, the pressure drop 
corresponding to a given constant flow rate rises strongly, 
leading to a sharp increase of the drag force F on the rod. 
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When the axis of the cyhnder is in the middle plane 
of the cell, the minimum hydraulic aperture of each of 
the two spaces between the cylinder and the nearest wall 
is ho{l — f3)/2 and the flows in both flow paths add up. 
If the rod is displaced from its equilibrium position and 
touches one of the walls along its full length, there is only 
one flow path of minimum hydraulic aperture /io(l — /3) 
(i.e. twice the previous one). In the lubrication limit 
(/3 1), the pressure drop Ap for a given flow rate q 
varies as the power —5/2 of the hydraulic aperture as 
shown by Eq. Even taking into account the fact that 
the flow takes place in a single channel, this variation 
with the aperture is fast enough so that the drag force is 
lower in the conflguration in which the cylinder touches 
the wall: more quantitatively, from Eq. ([3|), the coeffi- 
cient A_L decreases by a factor 2^/2 from its value in the 
centered position. This difference between the drag forces 
measured in these two locations of the cylinder is indeed 
observed experimentally as shown in Fig. [7] (dashed line). 

The case L ^ W seems therefore to correspond well 
to the 2D configuration studied numerically 0, [H, HO] ; 
however, when the cylinder is shorter than the width W 
of the slit, flow may be perturbed (particularly when /3 — > 
1) so that the 2D approximations do not reproduce well 
the observations. These effects will now be investigated. 



B. Variation of Ax with the length of the cyhnder. 

In this section, the variation of the transverse drag 
coefficient is studied as a function of the normalized 
length L/W characterizing the lateral confinement. 




FIG. 8: Variation of the normalized drag coefficient Ax/A^"^ 
as a function of the normalized length of the cylinder L/W , 
for cylinders of different diameters located in the middle plane 
of the cell, (a), (♦) and (•) symbols correspond respectively 
to experiments for /3 — 0.31, 0.65 and 0.83 (see Tab. HH for 
more details). Dashed-dotted and dashed line: 3-D numerical 
results obtained for /3 — 0.65 and 0.83. 

The experimental variations of Ax with W for three 
cylinders of different diameters d have first been com- 
pared. For each cylinder, the maximum value Ax"^ of 



A_L is reached when L/W is close to 1. In this case, the 
values of X™°-^ are close to the predictions of 2D simula- 
tions (see FiglT]) ; small fluctuations are observed and are 
likely due to experimental uncertainties (inhomogeneities 
of the cell aperture and cell roughness for instance). 

Figure [8] displays the experimental variation (symbols) 
of the normalized ratio A_l/A™"^ with L/W in the range 
0.42 < L/W < 0.99. As L/W decreases away from 1, 
the ratio Ax/A™"^ becomes significantly lower than 1; 
this variation occurs earlier and is particularly strong 
when (3 is large: for instance, for /3 = 0.83, A_l/A™°^ 
decreases by 60% for a small reduction of L/W by 5 %. 
This sharp variation is due to the partial diversion of the 
flow into the free space between the ends of the rod and 
the side walls which, in turn, reduces the drag force F. 
As observed in Fig.[8l this bypassing effect is particularly 
strong when the clearance at the end of the cylinder is 
large (low values oi L/W) while the interval between the 
surface of the cylinder and the two parallel walls is small 
(when /? ^ 1). 

In order to predict numerically the variations of 
X±/X±'^^, the full 3-D flow velocity fleld v(a;,y, z) must 
be determined, and not only the components Vy{y, z) and 
Vz{y,z) as before. This has been achieved by using the 
3-D version of the Freefem program (see Sec lIIip . As a 
validation test, we consider flrst the special situation in 
which the length of the cylinder is equal to the width of 
the cell (-L = W) . If, in addition, a perfect slip condition 
is used for the side walls, then the 3-D and 2-D simulations 
are equivalent and similar results should be obtained. Ac- 
tually, 3-D simulations predict values of Ax larger than 

2- D ones by 20 % or less: this difference is likely due to 
computational limitations related to the minimal practi- 
cal value of the mesh size. 

These 3-D simulations were then performed for different 
values of L/W and for two different normalized diame- 
ters /3 corresponding to actual experiments. The dotted 
and dashed-dotted lines in Fig [5] connect the data points 
corresponding to the ratios A^/A™"^ obtained from these 
simulations (the values of A™°^ are those obtained in the 
validation simulations for L — W). The experimental 
and numerical variations of A^/A™"^ with L/W are in 
very good agreement. This confirms that, in this geom- 
etry, the variations of the drag force with L/W reflect 

3- D modiflcations of the flow structure and cannot be ac- 
counted for by 2D models. 

The variation of Aj_ with the distance of the rod from 
the middle plane of the cell has also been investigated: as 
shown in Sec. IV Al the drag force should be lower when 
the rod is in contact with one of the front walls than 
when it is located halfway between them. We compared 
the experimental variations of Aj^/A™"^ with L/W in 
both configurations: the results obtained for two values 
of the normalized diameter /3 are displayed as symbols 
in Fig. [9] The drop of the coefficients when L decreases 
is more pronounced for a cylinder halfway between the 
walls. This can be explained by the different relative 
magnitude of the hydraulic impedance of the fiow paths 
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FIG. 9: Variation of the normalized drag coefficient Ax/A'"°^ 
as a function of the normaUzed length L/W for two cyhnders 
of different diameters located either in the middle plane of 
the cell (filled symbols) or near a wall (open symbols). (A, A) 
and (o, •) symbols correspond respectively to /3 = 0.31 and 
0.83 (see Tab. IIT|| . Lines are guides for eyes. 



between the ends of the rod and the side walls and of the 
direct paths between the rod and the front walls. 

Finally, the influence of the viscosity has been inves- 
tigated by comparing the results of experiments using 
identical parameters but fluids with two different viscosi- 
ties: rj — SOmPa.s (-I-) and r] = SOmPa.s (o). The points 
coincide, which confirms that the drag force is propor- 
tional to the viscosity. 



VI. CONCLUSION 

The experiments and numerical simulations reported 
in the present paper have allowed one to determine the 
influence of confinement effects on the drag force F on a 
static rigid cylinder in a viscous flow inside a rectangular 
slit. Significantly different results have been obtained in 
the cases of cylinders with their length parallel and per- 
pendicular to the mean flow, although, in both cases, F 
is proportional to the mean velocity and the viscosity in 
the range investigated and can, therefore, be character- 
ized by drag coefhcients A|| and A_l. 

In the parallel case, A|| increases linearly with the con- 
finement parameter /3 = d/liQ but does not diverge for 



/? = 1; increases much faster with /3 and diverges 
near (3—1 (due to the blockage of the flow) when the 
length L of the cylinder is close to the width W of the slit. 
2D numerical simulations in planes respectively perpen- 
dicular (for A||) and parallel (for Aj^ and L ~ W) to the 
mean flow reproduce well the results obtained in these 
two cases. In the perpendicular case analytical model 
based on the lubrication also provides a good agreement, 
still for L and for a strong enough confinement (i.e. 
for /3 > 0.2). 

When the cylinders are shorter than the width W of 
the slit, a bypass flow appears in the space between the 
between the edges of the cylinder and the side walls of 
the slit: this reduces the direct flow in the gap between 
the front walls and the cylinders. This effect is particu- 
larly strong when the confinement parameter (3 is close 
to 1; it results in a sharp decrease of the coefficient A^. 
3-D numerical simulations are needed in order to predict 
quantitatively this effect. 

The present experiments also provided evidence for an 
inertial lift force in the case of cylinders perpendicular to 
the flow direction. This force was observed for Reynolds 
numbers as low as 0.1 and kept the cylinders in the mid- 
dle plane of the model: this may explain recent obser- 
vation of the depinning of fibers trapped inside fractures 
[ssj . This observation may be contrasted with numerical 
simulations [U, HI] which predict an opposite effect, i.e. 
a force pushing the cylinder towards the walls. This lat- 
ter force should however appear only at larger Reynolds 
numbers: further studies are needed to investigate these 
issues. 

The present study dealt only with motionless rigid 
cylinders inside a viscous flow. Extending these studies to 
the measurement of forces on moving cylindrical objects, 
both rigid and flexible, will make the results applicable 
to the motion of freely swimming microorganisms [34] . 
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